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If Va or v T mix with v e , neutrino oscillations and collisions in a supernova (SN) 
core allow these flavors effectively to participate in (3 equilibrium and thus to obtain 
a large chemical potential. If a sterile species mixes with v e , these effects lead to an 
anomalous loss of energy and lepton number. We study flavor conversion in a SN core 
on the basis of a new kinetic equation which rigorously includes neutrino interference 
and degeneneracy effects. Our discussion serves as an example and illustration of the 
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properties of this "non-abelian Boltzmann equation" . 
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If neutrinos have masses and if the mass eigenstates do not coincide with the weak in- 
Q\ . teraction (or flavor) eigenstates, this "mixing" leads to oscillations and hence to transitions 
\ between different neutrino flavors: the individual flavor lepton numbers are not conserved. 
Q-V Therefore, in a thermalized neutrino ensemble, such as in a supernova (SN) core during and 
O ' immediately after collapse where neutrinos are trapped, all mixed flavors will be character- 
+3 • ized by the same chemical potential. Moreover, because the concentration of electron lepton 
c$ . number is initially rather large, the u e form a highly degenerate Fermi sea. The other flavors 
\ and v T are characterized by a thermal distribution at zero chemical potential unless they 
\ mix with v e in which case they would achieve the same large chemical potential. In a SN 
?h ' core heat and lepton number are transported mostly by neutrinos; the efficiency of these 
— ■ processes depends crucially on the degree of neutrino degeneracy for each flavor. Therefore, 
it is of great interest to determine the time it takes a non-z/ e flavor to equilibrate with v e 
under the assumption of mixing . 

If u e mixes with a sterile neutrino species, conversion into this inert state leads to the loss 
of energy and lepton number from the inner core of a SN. The observed SN 1987A neutrino 
signal may thus be used to constrain the allowed range of masses and mixing angles |§. 

Flavor conversion occurs by the combined action of two effects: Oscillations between 
neutrino flavors, and collisions with a "background medium" which break the coherence 
between the flavor components of a mixed state and thus interrupt the evolution of the 
phases. One may easily estimate the time scale for achieving chemical equilibrium to be 

t- 1 « sin 2 (26) J D (1.1) 
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where O is the neutrino mixing angle in the medium and D is an effective "interruption 
rate" for the coherent phase evolution, to be identified with a certain average of the weak 
scattering rates involving neutrinos @J5|]. Previous discussions [0,0,0 relied essentially on 
this reasoning. 

Eq. (|1.1|) is based on a single-particle wave function picture of neutrino oscillations and 
thus it is applicable if effects non-linear in the neutrino density matrices can be ignored. 
Therefore, this simple picture does not allow one to include rigorously the effect of Pauli 
blocking of neutrino phase space, an effect undoubtedly important in a SN core where the 
v e Fermi sea is highly degenerate. Therefore, we presently revisit the problem of neutrino 
flavor conversion on the basis of a newly derived kinetic equation || which rigorously includes 
neutrino degeneracy effects. Our discussion differs in several important aspects from previous 
works 03,0], but we agree with the overall picture and the numerical estimates to within 
factors of order unity. 

More importantly, our work serves as a detailed and explicit example for the use of the 
"non-abelian Boltzmann equation" that was derived in Ref. || and that we have extended to 
include neutrino absorption and production by the background medium. In this equation the 
oscillation aspects of neutrino propagation, which are usually treated as a 2-level Schrodinger 
problem, were unified with its kinetic aspects. Our present application is a particularly 
simple and transparent case because we may use the "weak damping" or "fast oscillation" 
approximation where the neutrinos can be treated as "propagation eigenstates" between 
collisions. We thus obtain a closed set of differential equations involving only the usual 
neutrino occupation numbers for the two mixing flavors. 

In Sect. II we proceed by casting the "non-abelian Boltzmann equation" of Ref. || into 
a form appropriate for the "weak damping limit". In Sect. Ill we study flavor conversion 
between v e and or u T , and we discuss the emission of sterile neutrinos from a SN core. 
In Sect. IV we summarize our findings. In Appendix A we discuss the special case of two 
mixing flavors with identical scattering amplitudes on a given target species and stress that 
flavor conversion does occur, contrary to statements found in the literature. In Appendix B 
we derive the neutrino scattering rates needed in Sect. III. 
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II. KINETIC EQUATION FOR MIXED NEUTRINOS 



A. General Result 



The usual objects of study for a kinetic description of a neutrino ensemble interacting 
with a thermal heat bath are the occupation numbers f p = (b Jf (p)b(p)) of a mode p of 
the neutrino field, where 6^(p) is a creation operator for a neutrino in this mode and the 
expectation value refers to the state of the ensemble, not to the vacuum. The physical 
situation of interest is one with highly degenerate neutrino Fermi seas so that we may 
ignore the presence of anti-neutrinos which otherwise would have to be included in a kinetic 
description of the system. 

The effect of neutrino mixing between n flavors can be rigorously included if we replace 
f p by a nxn "matrix of densities" Py(p) = (&j(p)&i(p)), i,j = l,...,n, where &|(p) is a 
creation operator for a neutrino of flavor i and momentum p The diagonal elements 

of p p = p(p) are the usual occupation numbers / p for each flavor i = 1, . . . ,n while the 
off-diagonal terms contain more subtle phase information. 

Earlier, Raffelt, Sigl, and Stodolsky || have derived a kinetic equation for the matrices 
p p of a homogeneous and isotropic ensemble under the assumption that all neutrino flavors % 
scatter on a given species of targets in the same way apart from an overall amplitude factor 
gi. If there are several species a of targets present, and if the corresponding quantum fields 
are uncorrelated, the collision integral is [notation dp = d 3 p/(2ir 
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pp = i\pp,n p ] d P' 



G a pp / G a (l- P p)+h. C . 
P pG a (l-pp,)G a + h.c. 

- yv PP - 2 



(2.1) 



where the nxn matrix G a = diag(g", . . . ,g°) in the neutrino interaction basis. Moreover, 
(<7") 2 Wp p / is the transition probability for a neutrino % with momentum p into a state p' 
due to the interactions with the medium species a. W a depends on the neutrino interactions 
with the medium constituents and on the intrinsic interactions of the medium. Finally, 
flp is a matrix of oscillation frequencies; its eigenvalues are the energies in the medium of 
the n neutrino flavors at momentum p. If there is only one flavor Eq. (|2.1|) reduces to the 
usual Boltzmann collision integral because p p — > f p , the oscillation term disappears, and G a 
becomes a simple number. 

Neutrino self-interactions are not included in this description. Because in a SN core the 
neutrino chemical potential is always much smaller than that of the electrons it is physically 
justified to neglect neutrino-neutrino interactions relative to other targets. 

Besides neutrino scattering we also need to include their absorption or production by the 
medium. Apart from pair-processes involving anti-neutrinos, which we continue to ignore, 
this can occur by charged-current reactions where a u e , for example, is transformed into an 
e~ or vice versa. We have derived the corresponding term for our non-abelian Boltzmann 
equation by following the same line of reasoning as in Ref. ||, except that we replace the 
neutral currents involving neutrino field bilinears with a charged current involving a neutrino 
and a charged-lepton field. 



3 



Our result can be represented in a compact form by means of the nxn matrices l{ which, 
in the interaction basis, contain only zeros except for a 1 on the i-th diagonal position. Then 
we find 

n 

/*„ = ... + £& - 1{4 P P }) v p - \{h, Pp } 4, (2.2) 
1=1 

where ". . ." stands for the terms already present in Eq. (|2.1| ) and { ■ , ■ } is an anti- 
commutator. As before, we assume that the effect of the neutrino interactions on the 
medium can be neglected so that V p , the production rate of a neutrino flavor i with momen- 
tum p, and «4p, the absorption rate of a neutrino % with momentum p, are given in terms of 
the equilibrium properties of the medium. For a single flavor the r.h.s. of Eq. (|2.2j ) takes on 
the familiar form of the difference between a rate-of-gain and -loss term, (1 — f p )V p — f p A p , 
including the usual Pauli blocking factor. 

We will always assume isotropy of the medium and of the neutrino ensemble so that all 
functions of p really depend only on E = |p|. Therefore, we define the angular average of 
the transition rate 

Wee> = ^ / Wp P '^ (2.3) 

where / dVL is an angular integration about the relative angle between p and p'. If we 
substitute this expression everywhere for W pp ' we may also substitute dp = <i 3 p/(27r) 3 — > 
E 2 dE/2ir 2 . Throughout most of our discussion we will prefer, however, to keep p and p' as 
our variables. 

Equations ( |2.1| ) and ( |2.2j ) are the "non-abelian Boltzmann collision integrals" which serve 
as a starting point for a kinetic treatment of neutrino flavor conversion in a SN core. 



B. Two-Flavor Mixing 

In what follows we shall focus on two-flavor mixing, n = 2. In this case the oscillation 
term in Eq. ( |2.1| ) can be simplified if we express all 2x2 matrices through Pauli matrices. 
Notably, we write p p = \{n p + P p • r) and Q p = E p + |B p • r, where n p = Tr(p p ) is the 
total number of neutrinos in mode p, E p = |Tr(f2 p ) is the average of the energy eigenvalues 
of this mode, and r is a three-vector of Pauli matrices. Oscillations alone then lead to the 
well-known precession formula 

P p = B p x Pp. (2.4) 

The "magnetic field" B p = B(p) around which the "polarization vector" precesses is given 
by the standard result 

£i(p) = l^sin2©o, 

^|P| 

/\ 771, ^ 

MP) = -^TT cos26) o - V2G F N, (2.5) 
2 P 
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where we have taken B 2 (p) = without loss of generality, Gp is the Fermi constant, Go is 
the vacuum mixing angle, Am 2 is the difference of the squared neutrino vacuum mass eigen- 
values, and we have used a relativistic approximation. Ignoring neutrino self-interactions 
which would provide only a small correction in the physical situation of interest we need to 
use N = N e — N^ T for the mixing between v e and For the mixing of f eiAt)T with a sterile 
species v x we have N = N e ^ T — \N n . 

The angle of the vector B against the 3-axis in flavor space is 20 p , twice the usual 
medium mixing angle, which is a function of the neutrino momentum. It is given by 

„^ sin20 o , n nS 

^ tan29 ^ CO3 2e - g / gr - (26 > 

where in our relativistic approximation E — |p| and 

Am 2 . . 

E r = 2.7 

2y/2G F N 

is a "resonance energy". For 9 < 1, a limit which we shall usually take, neutrinos with 
the energy E r experience maximum mixing with |0 P | = 7r/4. Besides t p we will also need 
s p = sin26 p and c p = cos26 p . We have 

S p = AQl + {l-E/E r f (2 ' 8) 

where, again, we have assumed O <C 1. 

We will mostly focus on the mixing of v e with some other flavor and will always assume 
that only z/ e 's can be produced or absorbed by the medium, while the other flavors can only 
scatter. Therefore, V p and A p without superscript will always refer to the production and 
absorption rate of u e . Of course, if a degenerate sea of v^s is produced by flavor conversion, 
the Fermi energy can exceed the muon production threshold, and in this case charged-current 
processes involving muons would have to be included. As an indication that we will exclude 
this possibility in our discussion we will use z/ e -z/ r -oscillations as our prime example. 



C. Weak Damping Limit 

Even for two-flavor mixing the general form of the collision integral remains rather com- 
plicated. However, for the conditions of a SN core we may apply a further approximation, 
the "weak damping limit" . It is easy to show that for the physical conditions of a SN core 
the oscillation rate for p p (the precession rate of P p ) is much faster than the scattering rate. 
Therefore, it is justified to consider density matrices p averaged over a period of oscillation. 
While the matrices p p are given by the four real parameters n p and P p which are functions 
of time, the matrices p p require only two such functions, for example n p and |P P |, while the 
direction of P p remains fixed and is identical with that of B p . 

One way of looking at the weak damping limit is that between collisions, rather than using 
weak interaction eigenstates, neutrinos are best described by "propagation eigenstates" , i.e., 
in a basis where the p are diagonal. In this basis the matrices of coupling constants G a are 
no longer diagonal, whence flavor conversion is understood as the result of "flavor- changing 
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neutral currents" . Of course, "flavor" now refers to the propagation eigenstates. However, 
because in general the effective mixing angle is a function of the neutrino energy one would 
have to use a different basis for each momentum state, an approach that complicates rather 
than simplifies the equations. Therefore, we will always work in the interaction basis. 

In order to derive an equation of motion for p we evaluate the r.h.s. of Eqs. ( |2.1|) and (|2.2|) 
in terms of the p p . We expand the result in Pauli matrices, leading to an expression of the 
form |(a p + A p • r). While in general the vector A p is not parallel to B p , the assumed 
fast oscillations average the perpendicular component of A p to zero so that only the parallel 
component remains. Therefore, 

j5 p = i[ap + (B p .Ap)(B p .T)], (2.9) 

where B p = B P /|B P | = (s p ,0,c p ). 

For our purposes the most practical quantities to parametrize the matrices p p are their 
diagonal entries, the occupation numbers of the two mixing flavors which we take to be v e 
and v T . It is straightforward to show that 

/p o ^ , (f;-f;K (o i 



+ ^f^(I o (2-10) 



where f^' T are the occupation numbers^ of 



For the charged leptons, we denote the occupation numbers by rep and n p . The chemical poten- 
tials and number densities are denoted by fj, UeT and N VeT for the neutrinos and fj, e>r and N e>r for 
the charged leptons. 
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From Eq. ( |2.9D we then find explicitly for the "absorption terms" Eq. (|2.2| ) of the collision 
integral 



J V 



Jp - 2 



2 ' /;, • ./„ x 



i 



i 



v. 



f e P + fl 



j p 1 j -p 



Ar 



(2.11) 



If neither v e nor i/ T are occupied because, for example, the medium is transparent to neutrinos 
so that they escape after production, we have = \s\V p , i.e., the production rate of v T is 
that of v e times | sin 2 (20), in agreement! with Ref. 0. 

In a SN core where the interacting neutrino species are trapped Eq. ( 2.11| ) is more 



complicated as backreaction and Pauli blocking effects must be included. It becomes simple, 
again, if s 2 1, because then the v e will reach f3 equilibrium much faster than the z/ r ; at 
this point we may use to lowest order the detailed-balance condition 



f e P {V p + A p )=V p . 



(2.12) 



Inserting this into Eq. (|2.11| ) leads to fl = |s 2 [(l — fp)V p — fpA p so that now the v T 
follow a naive Boltzmann collision equation with rates of gain and loss given by those of v e 
times | sin 2 (2B). If we consider the emission of a sterile species u x instead of u T , they would 
escape without building up so that their production rate would be \ sin 2 (26) that of u e , half 
as much compared to the above situation where the v e were also free to escape!. 



2 There is a small discrepancy with Turner's Q discussion who used an average production effi- 
ciency of sin 2 instead, which for small mixing angles is half of our result. He argued that, if 
v e was a superposition cos Qui + sin 1>2 of mass-eigenstates, the production amplitude of was 
sin© times that of v e and its production probability, therefore, sin 2 times that of v e . For small 
one has v T ~ V2, hence for small the production probability for v T appears to be sin 2 times 
that of v e . This line of argument, however, ignores the interference between v\ and V2 which leads 
to the oscillation phenomena. In the usual single-particle treatment the probability for finding a 
u T at a distance L from the production site of a v e is sin 2 (20) sin 2 (7rL/L osc ) where L osc is the 
oscillation length. Thus, on average one expects | sin 2 (20) neutrinos of the "wrong" flavor for 
each v e produced. 

3 In Ref. [Q] this factor 1/2 was not included, precisely compensating the missing factor of 2 
mentioned before. In Ref. 0] this factor was also left out. 
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For the scattering terms Eq. (p.l|) of the collision integral we find from Eq. (|2.9| ) explicitly 

fZ = \J dp'{ K'p/pMi - r P ) - >v PP '/;(i - [(4 - 4^ + (2 - 4)t P v^. 

+ [>Vp' P /^(l - £) - - fl*)] [s 2 p g 2 r - s p t P V^r 

+[>Wp'(i - /;) - w pp ,/;(i - /;,)] [4^ 2 - (2 - 4)v P '^r 
+[>Wp'(i - - vv#(i - #)1 fo e 2 + 4* P v^ 



(2-13) 



For notational simplicity we have written G a = diag(<?g, i.e., the indices e and r refer, 
again, to u e and z/ T . Moreover, we have suppressed the superscript a everywhere; a sum- 
mation over target species is understood. The equation of motion for f p is the same if we 
exchange e <-» r everywhere. In the absence of mixing we have s p = t p = 0, leading to the 
usual collision integral for each species separately. Integrating Eq. ( |2.13| ) over dp yields 



dp J dp' [>Vpp>(4 - 4) - f;,)(g 2 T + t p t pl g e g T ) 

-/p(l-/p')(^ + Vp'^ 

VW£(l-£)-W„,£(l-# 



(2.14) 



+ 



x 



s P 9e + s p'9 T ~ 2t p t p ,g e g T ) + (s p + s pl ) t p t p ,g e g T 



for the rate of change of the total v T number density due to elastic scattering. 

We will also have occasion to consider the mixing of v e with a sterile species v x for which 
g x = by definition. This yields for the scattering part of the collision integral 



(2.15) 



f; = f gl \ dp' w p/p f; f (2 -r p - r p ) - w pp ,(fZ + / P )(i - /pO 



If the v e stay approximately in thermal equilibrium, detailed balance allows for a further 
simplification. Together with the absorption term Eq. ( p. 11| ) and with Eq. ( |2.12j ) we find in 
this approximation 

r P = s \{ to - r P )v P - r P A p ] + g \ [ dp'\w p/p m - /*) - w pp ,fm - r P ,)] ) . (2.16) 



If the mixing angle is so small that the u x freely escape we may set f p = on the r.h.s. of 
this equation and find 



Nr 



dp^ 

y 4 



V p + I dp' g 2 e W plp r p , 



(2.17) 



for the emission rate of lepton number per unit volume. For the energy-loss rate we need to 
include a factor E — \p\ in the integral. 
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D. Small Mixing Angle 



In order to treat neutrino flavor conversion in a SN core it is not necessary to solve 
Eqs. (|2.11 ) and ( |2.13| ) in their full complexity. In a medium, maximum mixing obtains if the 
denominator of Eq. (|2.6| ) vanishes. For such conditions flavor equilibrium occurs practically 
on the time scale it takes to achieve (3 equilibrium which is instantaneous compared with 
other relevant time scales. In this case a detailed evaluation of the kinetic equations is 
superfluous. Therefore, we may focus on a situation where the neutrino mixing angle is 
small, 9 P <C 1. 

For small mixing angles (3 processes as well as elastic neutrino collisions with medium 
particles are much faster than the rate of flavor conversion which is suppressed by a factor B 2 . 
Therefore, both f£ and /p will be given to lowest order by Fermi-Dirac distributions which 
are characterized by slowly varying chemical potentials fJ>v T {t) and fi Ue (t). The integrated 
version of Eq. (|2.11|) with the detailed-balance condition Eq. ( |2.12|) , together with Eq. (|2.14|) 
then yields to lowest order in B 2 



dp&l\(i-f T P )v p -f;A p 



+ /dp/dp'Wpp' (e p 



glm-fl>)-glm-f\ 



+ (^B p -s e B P 2 /p(l-/ P r 
-(^e p -^B P 2 /p(l-/p e 0}, 



(2.18) 



and a similar equation for N Ue . Together with the condition of (3 equilibrium, ji n — fj, p — 
A*e — f^ve, that of charge neutrality, N p = N e , and the conservation of the trapped lepton 
number, d(N e + N Ue + N Vr )/dt = 0, these equations represent differential equations for the 
chemical potentials fi UT (t) and ii Va {t). 

Flavor conversion redistributes the lepton number trapped in a SN core among e~, u e , 
and v T . However, most of the lepton number resides in electrons and so, only a small change 
of the e~ and v e chemical potentials occur. Taking them at their unperturbed equilibrium 
values Eq. (|2.18|) simply is a differential equation for fi VT (t). 
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E. Very Degenerate Neutrinos 



In a SN core the v e Fermi sea is very degenerate. With regard to the neutrino distributions 
we may thus use the approximation T = so that neutrino occupation numbers are 1 
below their Fermi surface, and above. Also, elastic scatterings on medium particles will 
always lead to neutrino down-scattering because the medium can not transfer energy to the 
neutrinos: V^ee 1 — for E' > E. This implies immediately that the term proportional to 
(6p — 8p/) in Eq. (|2.18|) vanishes. Moreover, because /j, Ve > fi Ur the term proportional to 
/p(l — fp,) also vanishes. If we further observe that Eq. fl2.12| ) implies^ Ve = for E > fM Ve 
we find 



dE 



2tt 2 



+ dE'y.m, 



a 



'EE' 



f T e E - g a e e E ,) 2 E 2 E 

4vr 4 



12 



(2.19) 



where we have restored an explicit superscript a in the coupling constants and a summation 
over targets. Note that for degenerate neutrinos N u = /^/67r 2 whereas for electrons N e = 
fil/3n 2 because of the additional r.h. spin degree of freedom. Equation ( |2.19| ) is the starting 
point for our discussion of flavor conversion in a SN core. 



4 In Appendix B we will calculate Ve from the process e + p — > n + v e , taking the electrons to 
be highly degenerate, and the nucleons to be non-degenerate. In this limit one could be tempted 
to extend the integration up to fi e , i.e., to use Ve = for E > fi e rather than for E > fj, Ue , as 
was done, for example, in Ref. ||. This procedure, however, would violate the detailed-balance 
requirement Eq. ( |2.12| ) and thus, the approximations would not be self-consistent. In the final 
answer the upper limit of integration appears in a high power, the fifth power in the "vacuum 
limit" (see Sect. III). Therefore, even though fi Ue is not much smaller than // e in a SN core, /x^ e is 
much smaller than perhaps by as much as a factor of 10. 
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III. FLAVOR CONVERSION 



A. Mixing of v e with z/„ or v T 

We are now in a position to approach the problem of flavor conversion between the 
degenerate v e and or v T in a SN core quantitatively. To this end we consider two limiting 
cases for the small mixing angle. For a v e near its Fermi surface with E = fi Ue we have in 
Eq. O, taking N = N e , 

E _ (68keV) 2 1/3 4/3 

where P14 is the density in units of 10 14 gcm -3 and Yj gives the abundance of species j 
relative to baryons. If Am 2 is so large that E/E r <C 1 we may effectively use the vacuum 
mixing angle. In the opposite limit we may ignore cos2Bo in the denominator of Eq. ( |2.6|) . 
Hence we use 

icv i P "vacuum" 

l *l = 0oX U/£ "medium" ^ 

as our limiting cases. 

Given enough time, the v T will reach the same chemical potential as the v e . Therefore, 
it will be most practical to discuss the approach to chemical equilibrium in terms of a 
dimensionless v T density 

r] ^N l/ jN l/e . (3.3) 

Therefore, the equation of motion will be of the general form 

(97r) 1/3 G|iV e 5/3 F I(V) "vacuum" , 
fj = Q 2 x { (A m 2)2 a (3.4) 

1 ( f^I>^) "medium". 

For later convenience we have included the electron density N e rather than the baryon density 
as well as some numerical factors in the overall normalization. The F^' M are dimensionless 
functions of r]. It is interesting that in the "medium case" the dimensionfull factors do not 
involve Fermi's constant; a factor Gf? from the scattering rate cancels by Gp 2 from Q 2 E . 

Beginning with the first term in Eq. ( p,19| ) we use the neutrino production rate Ve due 
to the reaction p + e — > n + v e that was given in Eq. (Bl|). An explicit integration yields for 
the two limiting mixing angle cases defined in Eq. (|3.2|) 




(3.5) 



for the r.h.s. of Eq. ( |3.4j) . For comparison with Ref. we note that in the vacuum limit, 
and at t — when r\ = 0, our result can be written as N Ut = Oq (4G|/107r 3 )A^/i^ if we use 
C 2 + 3C 2 A = 4 
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Turning to elastic neutrino scatterings, i.e., to the term involving W EE i in Eq. (|2.19|) , 



we first consider nucleons as targets, v + iV — > iV + v. This neutral-current (NC) reaction 
involves identical amplitudes for all flavors whence our matrices G can be taken to be the 
2x2 unit matrices for both proton and neutron targets. Then Eq. ( |2.19| ) implies immediately 
that no flavor conversion occurs in the limit of vacuum mixing where &e — @£" = ©o! hence 
Fnc{v) = 0- m the limit of medium mixing where oc E~ l , however, NC reactions do 
contribute to flavor conversion, contrary to statements found in the literature |IJ, an issue 
discussed in Appendix A. By means of the function given in Eq. (TBI) we find explicitly 
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]=p,n 



where A p n are the density fractions of protons and neutrons with X p + X n = 1. The flavor 
conversion rate is much smaller than that of CC scattering Eq. (|3.5|) because it is suppressed 
by the "recoil factor" (fi Ue /m N ) 2 . Put another way, damping is small because the neutrino 
energy in a collision does not change very much so that even in the medium Qe ~ Qe 1 - It is 
nevertheless conceptually interesting that this process alone would suffice to reach chemical 
equilibrium. 

Next, we consider r.h. electrons as targets. According to Eq. (|BTJ) we have gf" = ~ ~ 
so that no flavor conversion occurs in the vacuum case: F^(t]) = 0. For the medium case 
we find with the approximate transition rates of Eq. ( |B12|) 



F * {r]) = 1 7? (® ~ ^ 2/3 + 77 + > giv) ~ ' ]r/4/3 ~ ^ 5/3 )- (3 - 7) 

Chemical equilibrium corresponds to t] — 1; it is easy to check that F^(l) = 0. 

For l.h. electrons we have according to Eq. (p6|) g\ ~ 3/2 and g\ ~ — \. Therefore, we 
need to consider the vacuum as well as the medium case. We find with Eq. ( |B12|) 

/•£■(„) =rA — } (1 - ^ 1/3 ) 3 (l + b 1 " + ToV 2/3 + h 




F^(V) = ^ ^ (lv- 1/3 + \ logfa) - 1 + IV 1/3 + h 2/3 - Sj^ 4/3 ) • (3- 



32 fi 

In the medium case the r.h.s. of the equation of motion diverges initially when r] = 0. This 
singularity is unphysical because it is related to the mixing angle oc E^ 1 becoming 
"infinite" at low neutrino energies; in reality \Qe\ never exceeds 7r/4. 

In the "medium case" the contribution of r.h. electrons is much smaller than that of l.h. 
ones, in part because the function Wee' is much smaller for r.h. states (see Fig. ||), and in 
part because for r.h. electrons the integrand (g T ®E — ge®E') 2 oc (E — E') 2 /E 2 E' 2 is much 
smaller than for l.h. states where we have (g T QE ~ 9e®E') 2 oc (3E + E') 2 / E 2 E' 2 . In any 
event, we may now add the two contribution and obtain 

FrM+l(7/) = h (^~ 1/3 - W + K 1 + ^ 4/3 ) Mu) + h 1/3 + v - & 4/3 - ( 3 - 9 ) 

for the total electron contribution. 
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We now collect our results. In the vacuum limit the only contributions are from CC 
nucleon scattering, see Eq. (|3.5|), and from the scattering on l.h. electrons, see Eq. 



The electron contribution is suppressed by a relative factor ^(/^//^e), and it has a different 
structure of the r\ dependent term. In Fig. [I] (upper panel) we show these terms; the electron 
contribution drops to zero much faster than the nucleon term for rj — > 1. This is due to 
degeneracy effects because as the v T Fermi sea fills up the accessible energy difference E — E' 
between initial- and final-state neutrinos decreases so that an ever smaller fraction of the 
electrons is available as scattering targets. For all practical purposes we may neglect the 
electron term entirely so that the deviation from flavor equilibrium relaxes as 

|(l-^) = - iL ^ Lj (3-10) 
at r 

where 

In the medium limit we have a contribution to the damping rate from CC scattering on 
nucleons, given by Eq. ( |3.5| ). NC scattering also contributes, see Eq. ( |3.6| ), but this term is 
suppressed by a small "recoil factor" so that we may ignore it. NC and CC current scattering 
on electrons both contribute; the combined result was given in Eq. ( |3.9|) . However, compared 
with the nucleon result it is suppressed by a factor ^(/v/aO, while the r] dependent terms 
are shown in Fig. [L] (lower panel). Again, the electron contribution quickly drops to zero 
as flavor equilibrium is approached. Thus, while electrons initially contribute to flavor 
conversion, they may be neglected for the overall time scale and we may use, 

^(l-,) = -(l^) (3.12) 



where 



2vriV e , / IkeV \ Y e p 



T Qf.= — * = 1.8xlO~ J s ; - A —. — -. (3.13) 

(Am 2 ) 2 \y/A^J 10 14 g/cm 3 v ; 

In Fig. ^] we show contours of log 10 (rOo) according to these results. Considering the 
allowed range of vacuum mixing angles it is evident that neutrino mass differences well in 
excess of the cosmological bound of around 100 eV are required to achieve flavor equilibrium 
on a time scale of seconds or faster. 



B. Sterile Neutrinos 

As another application of our equations we may easily calculate the emission of "sterile 
neutrinos" v x from a SN core if these particles mix with u e , a problem previously discussed 
in Ref . . The v x would escape directly from the inner SN core so that the relevant figure of 
merit is the SN energy-loss rate Q, and the loss rate of lepton number Nl into this channel. 
A general expression for Nl was given in Eq. ( |2.17| ). 
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We consider, again, the two limiting cases for small mixing angles defined in Eq. ( |3.2|) 
and begin with the "vacuum limit" where Am 2 is so large that typically E <C E r . Assuming 
a small vacuum mixing angle we may thus use s 2 = 4Go in Eq. ( |2.17|) . We first consider the 
CC process e + p — > n + v and find with Eq. (|B1|) 



Nl= e 2 ^G 2 F N ef t 5 



1 



Q = 6 2 — G 2 iV e/ < (3.14) 

where we have used Cy + "5C\ ~ 4. Similarly, we consider the NC scattering process 
v + iV — > N + v which now does not involve any destructive interference because g x = 
for all processes. Therefore, we obtain the same result as in Eq. ( |3.14j ) if we incorporate 
an extra factor (N p + N n )/N e = 1/Y e , and a further factor | for the reduced coupling 
strength of NCs relative to CCs. (For NC scattering we also take C v + "iC\ « 4 for both 
protons and neutrons.) For electrons we use the approximate transition probabilities given 
in Eq. ( |B12|) and add the contributions from r.h. and l.h. electrons from the start, recalling 
that {g^) 2 = j while (g,P) 2 = §. Then we find the same results as Eq. ( |3.14j ) except for 



an extra factor ^(fi Ve / fi e ) for Nl and ^(^i Ue / fj, e ) for Q. Clearly, we may neglect the 
contribution of electrons relative to those of nucleons. Therefore, lepton number is lost at a 
rate 

Y L = -GlGlNF 12(36 5 ?r)1/3 (Y e + \)Y u f. (3.15) 
Because Yl ~ Y e and Y Ue is Y e up to a factor of order unity, we have 

^3 « -elGlNf = -e 2 a 2.7X10'" jp^. (3-16) 

Since most of the trapped energy is stored in the degenerate leptons, energy is lost on a 
similar time scale. In order to avoid conflict with standard SN physics, and especially with 
the observed duration of the neutrino signal from SN 1987A, we must demand 0o ^ 10 -5 , 
a bound in agreement with that found in Ref. ||. 

Next, we turn to the "medium limit" defined in Eq. ( |3.2|) where Am 2 is small enough 
that typically E ^> E r . For degenerate v e this means fx e ^> E r or (Am 2 ) 1 / 2 < 100 keV for the 
inner core of a SN. It is easy to evaluate Q and Nl f° r the CC and NC processes involving 
nucleons. However, in the present situation the scattering term on electrons dominates 
because Q and Nl diverge if we use the simple approximation 6^ = Q E r /E of Eq. (ft.2|) . 
We encountered this divergence before in Eq. ( p.8|) , but there we were able to ignore it 
because it affected only the initial process of flavor conversion while the overall time scale 
was dominated by the CC nucleon scattering process. In the present case, the v x Fermi sea 
never fills up because they continuously escape so that the resonance region with E' = E r 
always remains within the region of integration. 

In order to obtain a meaningful result we need to keep the full expression Q2.8| ) for s 2 . 
However, this means that the range of integration includes a certain regime where the mixing 
angle in the medium is so large that even the "sterile" neutrinos are trapped. There will be 
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a critical mixing angle C below which the sterile neutrinos escape freely, and above which 
they are essentially trapped. We assume that 9q C 6 C < 1 so that finally 



where E c = E r (l + 9 O /0 C )- An explicit integration yields, keeping only the lowest powers 
of 6 and of E r /fi Ue , 



9 GlN e al E r 

N L = 9 9 c -— r . 3.18 

807T d \l e 



Because of the resonance, the emission rate is proportional to 0oO c rather than to Oq. 
The energy- loss rate is found by including an extra factor of E' in the integral Eq. (|3.17|) . 



To lowest order we find Q = N^E r because of the resonance. Put another way, in the present 
situation mostly neutrinos are emitted with energies around E r <C fi Ue . Hence lepton number 
is lost much faster than energy. After the lepton number of the SN has been depleted by this 
mechanism, assuming it dominates the standard diffusion transport, most of the energy will 
still be trapped and has to be emitted by diffusion to the surface. Therefore, a constraint on 
Am 2 and O would require a much more detailed consideration of the cooling history of a 
SN core with an early loss of lepton number. In this regard we disagree with the conclusions 
of Ref. |J where even in this case a bound was derived on the basis of a simple energy-loss 
argument. There, the scattering on electrons and thus the quick deleptonization was not 
included. 
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IV. CONCLUSIONS 



Starting with the general and rigorous "non-abelian collision integral" for mixed neutri- 
nos of Ref. we have derived a simple and transparent kinetic equation for the occupation 
numbers of two mixed neutrino flavors under the assumption of fast oscillations relative to 
the collision rate, an approximation appropriate for a SN core ("weak damping limit"). This 
form of the equation is both suitable for practical calculations concerning neutrino flavor 
conversion in a SN core, and also allows one to develop a better understanding of the inter- 
play between collisions and oscillations encapsuled in the relatively complicated expression 
of Ref. §. 

We explicitly derived a number of simple general results, some of which are in conflict 
with statements found in the literature. It was thought, for example, that two neutrino 
flavors could not reach chemical equilibrium with each other if for each target species they 
had identical scattering amplitudes. In this situation it was supposed that the medium could 
not "measure" the flavor composition of a mixed state. We found that this conclusion is 
correct only if the mixing angle is independent of the neutrino energy. As another example, 
we found that the emission rate of sterile neutrinos which mix with v e with an angle are 
emitted with a rate | sin 2 (20) times the production rate of z/ e 's if the medium is transparent 
to the z/ e 's. On the other hand, if the u^s are approximately in thermal equilibrium, the 
relative rate of gain for the sterile states is only \ sin 2 (20). 

We used our collision integral to calculate the flavor relaxation time r between the 
degenerate v e and the non-degenerate or v T in a SN core. A contour plot of tOq (vacuum 
mixing angle ©o) was shown in Fig. ^| as a function of the neutrino squared mass difference 
Am 2 , and of the density. Neutrinos obeying the cosmological mass bound can not achieve 
chemical equilibrium on a time scale below that of SN cooling. 

We also calculated the time scale for the emission of the stored energy of a SN core if v e 
mixes with a sterile species v x . If VAm 2 > 100 keV the vacuum mixing angle was found to 
be constrained by 0o ^ 10~ 5 . If V Am 2 <C 100 keV the "resonance energy" where u e and v x 
become degenerate is below the v e Fermi surface, leading to quick deleptonization by virtue 
of the scattering process v e + e — > e+v x where the final-state v x carries the resonance energy. 
In this case, no simple constraint on Bo can be derived. 
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APPENDIX A: NEUTRAL-CURRENT SCATTERING 



There is a limiting case of our kinetic equation that deserves some special attention. 
Neutral-current (NC) neutrino interactions cause a large contribution to the scattering rate 
in a SN core. Do they also contribute to flavor conversion? In order to answer this question 
we first note that the relevant G a matrix in Eq. ( |2.1| ) is the unit matrix because the NC 
scattering amplitude is the same for all neutrino flavors; the corresponding transition rate 
is denoted by W^. Therefore, the collision integral is 

Pp = i[p P , flp] + / dp'[w™( Pp , - \ {/v, Pp }) - < p S( Pp - §W, Pp })] , (Al) 

where we have ignored all other processes. For the global density matrix p = J p p dp we 
thus find 

p = i Jdp[p p ,Q p ] (A2) 

while the collision term vanishes identically. Therefore, it appears that NCs do not contribute 
to flavor conversion 

Of course, some flavor conversion takes place by virtue of the oscillation term Eq. (|A2|) . 
As a first case we assume that in the medium the mixing angle is dominated by the vacuum 
value Bo which is the same for all modes p. Then the precession of the polarization vector 
of all modes is around the same direction. If we start with a pure ensemble of v e the final 
v T population is |sin 2 (26 ) of the total neutrino density. This result remains valid even if 
the neutrinos are scattered around between different modes: there is no flavor conversion 
beyond that caused by free oscillations. 

This situation is better understood in the weak damping limit where oscillations are 
fast compared with the rate of collisions. Then, the trivial free oscillation effect takes place 
"instantaneously" after preparing the ensemble. The subsequent change of the total v T 
density is found by setting g e = g T = 1 in Eq. Q2.14| ), 



K = Jd P J d V ' w^{(4 - Sp ,)(i + t P v)[/p r (i - r pl ) - /p(i - m 

+K + 4 - t P v(4 + 4o][#(i - %) - r p {\ - #)]}. (as) 



It is evident that Eq. ( |A3[ ) vanishes identically if the mixing angle does not depend on p. 

This picture changes dramatically, however, when the medium mixing angle is a function 
of p. In this case N Ue and N Ut do not, in general, vanish. This behavior is understood if we 
recall that in the weak damping limit the neutrino states between collisions are most natu- 
rally described as "propagation eigenstates" . As both flavor components scatter identically 
by NCs one arrives at the naive conclusion that collisions cause no flavor conversion because 
the coherence between the two components is not broken. However, if the momentum state 
changes in a collision, the mixing angle of the new state is different from the original one 
so that the new propagation eigenstate is a different combination of flavor eigenstates. The 
subsequent adjustment by the assumed fast oscillations leads to a certain amount of flavor 
conversion in each collision. 

This behavior can be made more explicit in the limit of small mixing angles where 
Eq. ( |A3| ) becomes 
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K = Jd P J d P 'w^{(e 2 p - e p ,)[/ p r (i - %) - r v {i - #)] 

+(©p - ©p0 2 [^(i - /;) - /p r (i - m}- (A4) 

Moreover, in the limit 0-^0 the distributions / p and / p are in kinetic equilibrium so that 
with a non-zero but small G 2 they will deviate from equilibrium only to this order. Hence, 
we may use the detailed-balance condition 

/ dp' [W p N p S/ p (l - / P - W^f p/ (l - / P )] = 0(© 2 ) (A5) 

for both v T and v e . Next, we multiply this equation with 2 , apply / dp, and substitute 
p «-> p' in the second term under the integrals. Therefore, the first line in Eq. (|A4[) is of 
order 9 4 and may thus be neglected: 

n Ut = J dp J dp' w p N p S(e p - e P o 2 [^(i - #) - f£(i - #)]. (A6) 

If the ensemble consisted originally of v e only, and if the distributions are in kinetic equilib- 
rium to order O 2 , then / p > fl + (9(0 2 ) until chemical equilibrium is reached. Hence the 
r.h.s. of Eq. ( |A6| ) is always positive until the process of flavor conversion is complete. 

Therefore, contrary to the naive conclusion stated above, flavor equilibrium will be 
achieved by NC scatterings alone as long as 6 P 7^ 6 p / for |p| 7^ |p'|. This result, too, 
must be interpreted with care, however. If NC interactions were the only interactions there 
would be no difference between the refractive indices of the two flavors. Then, the mixing 
angle in the medium would be strictly identical with that in vacuum and thus independent 
of p; hence flavor equilibrium can not be achieved. 

Thus, our result really means that CC interactions have a two-fold effect on the process 
of flavor conversion. They cause different scattering amplitudes for the two mixing flavors 
and thereby break the coherence between the components of a mixed state, leading to 
decoherence. Also, they cause different forward scattering amplitudes and hence lead to a 
medium-induced component of the oscillation term. This implies that the mixing angle in the 
medium becomes a function of the neutrino energy, allowing NC interactions to contribute 
to flavor conversion. 

Of course, in a SN core the dominant NC scattering rate is on nucleons which, because 
of their relatively large mass, do not allow for a large amount of energy transfer in a collision 
with a neutrino. In this case the neutrino energy before and after a collision remains ap- 
proximately the same so that P ~ p / . In this case the impact of NCs on flavor conversion 
is relatively small as we saw quantitatively in Sect. III. 



APPENDIX B: NEUTRINO INTERACTION RATES 

In order to evaluate the equations of motion for the neutrino distributions in a SN core 
we need the production rate Ve of electron neutrinos with energy E, and the transition 
probability W EEI for the scattering on a target species a of a neutrino with energy E into a 
state with energy E'. 
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1. Charged-Current Scatterings on Nucleons 



Beginning with Ve, electron neutrinos can be produced and absorbed by the charged- 
current (CC) reaction p + e <-> n + v e . Because the effect of flavor conversion causes only 
a small perturbation of the electron distribution we compute Ve from the medium equi- 
librium properties. Moreover, we work in the limit of highly degenerate leptons so that 
the occupation numbers of the electrons are taken to be 1 below their Fermi surface, and 
above. The nucleons, however, are only partially degenerate; as a first approximation 
we ignore nucleon Pauli blocking factors, and we ignore correlation effects in the medium. 
Finally, the nucleons are non-relativistic so that the CC scattering of an electron of energy 
E produces a neutrino of the same energy. Therefore, we approximate Ve = &eN p where 
N p is the proton density and oe is the CC scattering cross section for electrons of energy 
E; it is o E = (C v + 3C^)GpE 2 /tt where Cy and Ca are the usual vector and axial weak 
couplings which in vacuum are 1 and approximately 1.26, respectively. The exact value of 
Ca in a nuclear medium is not known, but likely it is suppressed somewhat so that 1.0 is a 
more realistic value ||. Therefore, we use (Cy + 3C|) ~ 4 and find altogether 

V E *^E 2 N p , (Bl) 

7T 

accurate up to a factor of order unity. (Pauli blocking of nucleons can not be neglected 
entirely; a suppression by up to a factor of 3 or so occurs.) 



2. Neutral-Current Scatterings on Nucleons 



Turning to elastic neutrino scatterings we first consider nucleons as targets, u + N 
N + v. Because the nucleons are non-relativistic we may use 

27r 2 da Pt1l 



"'EE' 



<P,n E , 2 dE , (B2) 

where a p ^ n is the scattering cross section on p or n for neutrinos of energy E. A standard 
calculation yields M 



da 
dE' 



Glm N 
2tt 



A + B i-\ +C^^m N 



(B3) 



where A = \{Ca + Cy) 2 , B = \{Ca — Cy) 2 , and C = \(C\ — Cy). In vacuum, for protons 
Cy = 1 — 4sin 2 Bw and Ca = 1-26 while for neutrons Cy = —1 and Ca = —1.26. In a 
nuclear medium the axial charges for NCs are likely suppressed by a similar amount as for 
CCsi. 

Because the nucleons are relatively heavy compared with the neutrino energies we go 
only to lowest order in E/m^ where is the nucleon mass. In this approximation we find 



VV EE> 



4E 2 



j=p,n 



2(C 2 A + C 2 ) 3 + (C 2 A - C 



E-E' 



v)j 



E 2 



m N 



(B4) 



To this order the maximum energy transfer in a collision is 2E 2 /m^ so that Wee' — for 
E' > E and for E' < E(l - 2E/m N ). 
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3. Neutrino Electron Scattering 



The scattering of v e and v T on electrons can each be described by an effective NC Hamil- 
tonian. Taking the neutrino field i\) v to represent a two-spinor in flavor space we have 



H^ R = £| ^(1 T 7 5 )^e VV7 M (1 - 7 5 )G L 'V- (B5) 

The upper sign refers to L (left-handed electrons), the lower to R (right-handed electrons), 
and 

L= /2sin 2 e w + l 
' \ 2 sin 2 6 W - 

GR= /2sin 2 6 w \ 
V 2sm 2 e w / 

where Bw is the weak mixing angle. We will always use the approximate value sin 2 0\y « 
1/4. 

Because the electrons are ultra-relativistic the chirality states coincide essentially with 
the left- and right-handed helicity states. Therefore, the two electron helicities appear as 
two separate target species. Treating the electrons as massless free Dirac particles we find 

W - = /l7/^r/(2^/(2^ 2g2g-2^- " h(1 " "*> W ^ + k -" ~ k) (B7) 

where p and p' refer to the neutrinos, k and fc' to the electrons, and are the electron 
occupation numbers. The squared matrix elements are |.M L | 2 = 8Gp(pk)(p'k') and |.M R | 2 = 
8Gp(pk')(p'k), respectively ||10|| . Performing all angular integrals explicitly we obtain 

(7 2 /?3 r oo 

Wee' = x£e*J q dujn ^ 1 ~ n ^ A ^I E )^ 

r<2 77/3 „oo 

W * E ' = 30^ Wo dUJnM ~ U ^ AR{UJ/E ^ (B8) 

where is the electron occupation number at energy uj (we have assumed isotropy of the 
medium) and uj' = uj + E — E' . Moreover, we have 

A ( °-\ioe 2 ±5e + i if e > i, ( 9) 

where the upper sign refers to L, the lower sign to R. 

We stress that Wee' is different for l.h. and r.h. electron targets because the angular 
momentum budget of these reactions is different. In the CM system v-e^ scattering is 
isotropic, while v-e R scattering has a (1 — cos^) 2 angular distribution, i.e., it is forward 
peaked ||. This implies that in the system of the background medium v-e scattering is 
more efficient at transferring energy to l.h. electrons than to r.h. ones. 

In order to perform the uj- integration we take the electrons to be highly degenerate, i.e., 
we take = 1 for uj < fi e and otherwise. Moreover, we take the neutrino energies to be 
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less than /i e , and we note that neutrinos can only down-scatter; hence we have \i e > E > E'. 
Thus, only electrons with energies above \i e — E + E' contribute and we have 



du n^il — n^i) . . . — > / duo ... . (BIO) 

J^-E+E' 



Then we find 



W%b> = Wo (E E E 2 ' W W\Ej^ E'l^\ 

{E ~ E 2 ' )E ' W^E/^E'/^) i£2E<» e + E', 
W\e> = W x\ / (Bll) 

W 2 L {E/fi e , E'/n e ) if 2E > fi e + E', 



E 2 E' 2 

where W = nGpN e /fi e and 

W R (x, y) = i[60 - 30(2x - y) + (20x 2 - 25xy + lly 2 )], 
Wfo,y) = i[60 - 30(a; - 2y) + (llx 2 - 25xy + 20y 2 )}, 
W^(x, y) = - 6y - 15y 2 + 20(2x 3 - y 3 ) - 15(2x 4 - 4x 3 y + y 4 ) 

+6(7x 5 - 15x 4 ?/ + 10xV - y 5 ) 

-29s 6 + 36a; 5 ?/ - A5x 4 y 2 + 20x 3 y 3 - y 6 ]. 

These expressions are already fairly complicated, and the number of terms multiplies 
enormously once we perform the integrals over dE' and dE to obtain the damping rate. 
Therefore, we include only the leading terms in Ej\i e and E'/fi e since the neutrino energies 
are always smaller than the electron chemical potential. To this order we use 

w r _ w Cg - E')E' 
Wee 1 — w o ^ , 

Wee> = W {E ~ E E ' )E . (B12) 

In Fig. |] (upper panel) we show (W^ E ,/Wo)(E' / E) 2 for E = i/x e , |/x e , |/x e , and /x e as a 
function of E'/fi e . The solid lines refer to the exact result Eq. ( Bll ) while the dotted lines 
refer to the lowest-order result Eq. ( |B12|) . In the lower panel we show the same graphs for 
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FIGURES 



FIG. 1. The 77-dependent part of the functions F^(r)) (upper panel) and F^(rj) (lower panel). 
The curves from CC scattering on nucleons are according to Eq. (|3.5| ), while for electron scattering 
the "vacuum" curve is from Eq. (|3.8|), the "medium" curve from Eq. (p.9|). 



FIG. 2. Contour plot for t0q according to Eqs. ( |3.11|) and ( 3.13 ), taking /i„ e / ' \i e m 1. The 
contours are marked with log lo (r0g). 

FIG. 3. Transition probability Wee' for i^-e-scattering for r.h. electrons (upper panel) and l.h. 



electrons (lower panel) . Solid lines: Exact result according to Eq. ( Bll ). Dotted lines: Lowest-order 
approximation according to Eq. ( |B12| ). 
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